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While quantum mechanics exquisitely describes 
the behavior of microscopic systems, one ongoing 
challenge is to explore its applicability to systems 
of larger size and mass. Unfortunately, quan¬ 
tum states of increasingly macroscopic objects 
are more easily corrupted by unintentional mea¬ 
surements from the classical environment. Addi¬ 
tionally, even the intentional measurements from 
the observer can further perturb the system^. 
In optomechanics^, coherent light fields serve as 
the intermediary between the fragile mechani¬ 
cal states and our inherently classical world by 
exerting radiation pressure forces and extract¬ 
ing mechanical information. Here we engineer 
a microwave cavity optomechanical system® to 
stabilize a nonclassical steady-state of motion 
while independently, continuously, and nonde- 
structively monitoring it. By coupling the mo¬ 
tion of an aluminum membrane to two microwave 
cavities, we separately prepare and measure a 
squeezed state of motion^. We demonstrate a 
quantum nondemolition (QND) measurement®"^ 
of sub-vacuum mechanical quadrature fluctua¬ 
tions. The techniques developed here have direct 
applications® in the areas of quantum-enhanced 
sensing® and quantum information processing, 
and could be further extended to more complex 
quantum states^®. 

Reflecting light off a mechanical object induces a mo¬ 
mentum transfer, allowing one to control and measure 
the mechanical state. When the photon scattering rate 
exceeds the phonon decoherence rate the mechanical sys¬ 
tem becomes more strongly coupled to the photon reser¬ 
voir than to its own thermal environment. This regime 
is usually obtained by embedding a mechanical resonator 
into an electromagnetic cavity to increase the interac¬ 
tion strength per photon^. Additionally, the cavity fil¬ 
ters the density of states available for the scattered pho¬ 
tons, allowing control over the ratio of Stokes and anti- 
Stokes scattering rates. Importantly, the nature of the 
optomechanical interaction implies that the light held 
interacts with both mechanical quadratures, with fun¬ 
damental consequences on the mechanical state prepa¬ 
ration and measurement. On one hand, the precision 
on the simultaneous measurement of both mechanical 
quadratures is limited by the Heisenberg’s uncertainty 
principle^^fl^. On the other hand, the state preparation 
via sideband cooling exploits the coherent exchange of 
the cavity and mechanical state and is therefore limited 
by the statistics of the classical light field^®"^®. 

Both limitations can be overcome using polychromatic 


coherent light. One can address and manipulate each 
mechanical quadrature differently by engineering inter¬ 
ference processes between their couplings to the cavity 
quadratures. More specifically, two-drive schemes can be 
used to design a single quadrature measurement of the 
mechanical oscillator, known as backaction evading^^"^. 
The scheme fulfills the requirement for a QND measure¬ 
ment, which is an important tool for the tomographic 
reconstruction of arbitrary quantum states. A similar 
scheme was proposed by Kronwald at al'^ to prepare a 
mechanical squeezed state, following an analogous idea 
formulated for trapped ions^®’^^. In this work, we si¬ 
multaneously implement these two schemes in a single 
microwave optomechanical system to perform the tomo¬ 
graphic measurement of a squeezed state of a macroscopic 
mechanical oscillator. 

Consider a cavity optomechanical system where the 
position of a mechanical resonator of frequency Clm tunes 
the resonance frequency lOc of an electromagnetic cavity^. 
Two drives are applied to the cavity, at both mechanical 
sidebands = Wc ± The strength of each drive can 
be parametrized by its scattering rate T^ = j k, 

where go is the vacuum optomechanical coupling rate, k 
is the cavity linewidth and is the number of intra¬ 
cavity photons induced by each drive. Assuming a me¬ 
chanical relaxation rate F^ and the condition F^, F^ ^ 
K <C Ulrm one can write the relations between the mechan¬ 
ical quadrature amplitudes, Xi and X 2 , and the ampli¬ 
tude and phase quadrature of the cavity fields, A and ip, 
reading"^ : 


(Xf) 

in) 
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Here {X^^ = 2n^ -I- 1 is the variance of the mechan¬ 
ical quadratures for an equilibrium thermal occupancy 
n^, and {p^) = (A^) = 1 are the variances of the cavity 
quadratures for an ideal coherent state. For F+ = 0, cor¬ 
responding to driving only the lower sideband, one recov¬ 
ers the sideband cooling limit and at high scattering rate, 
F" Tm, each quadrature of the mechanics is cooled to 
the cavity quadratures, (Af) = {gA) and (A|) = (A^). 
Another limit is F"'" = F", corresponding to driving sym¬ 
metrically the upper and lower sidebands. This is the 
case of a drive on resonance with the cavity whose am¬ 
plitude is modulated at a mechanical frequency, perform¬ 
ing a QND measurement of the mechanical quadrature 
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FIG. 1: Device description and experimental setup, a, False-colour optical micrograph of the aluminum device (in grey) 
on a sapphire substrate (blue). Centered at the bottom of the micrograph is a mechanically compliant vacuum gap capacitor. 
The capacitor’s electrode is split into two plates, each shunted by a different coil inductor, giving rise to two microwave LC 
resonances, b. Frequency space diagram. Except for the mechanical linewidth, all the frequencies and linewidths are to scale. 
The microwave cavities have Lorentzian densities of states of width Ki/27r = 1.7 MHz and K2/27r = 2.1 MHz, centered at 
a;i/27r = 8.89 GHz and a;2/27r = 9.93 GHz. These resonance frequencies are tuned by the motion of the top plate of the 
capacitor, at the mechanical frequency flm/27r = 14.98 MHz. The red and blue dashed lines indicate the four mechanical 
sideband frequencies at which the cavities are driven, uj^2 = ^1,2 ± c. The circuit is placed on the cold stage of a cryogenic 
refrigerator (base temperature T = 30 mK). Up to four strong microwave drives and one weak microwave probe are inductively 
coupled to both cavities via a single port. The reflected signals and the noise emitted by the cavities are cryogenically amplified 
and demodulated at room temperature. 


Xj. Indeed, under these conditions, Eqs.l and 2 read 
(1?) = (X^2J (X|) = (X,l) + (4r-/r^)(i2). The 

Xi quadrature is unaffected by the measurement and 
the radiation pressure shot noise backaction is placed on 
the orthogonal quadrature X 2 . Finally the preparation 
of a squeezed state occurs in the intermediate regime, 
r+ < r . The mechanical mode is coupled, at a reduced 
rate T” — to an effectively squeezed microwave bath, 
whose minimum variance is (\/r“ — •\/r+ )^/(r~ — T"*') < 
1 . 

In order to separately prepare and read out a me¬ 
chanical state, we engineer a microwave optomechani¬ 
cal system where a single mechanical mode is coupled 
to two microwave cavities. The experimental setup is 
shown in Fig.I. The circuit, made out of aluminum on 
a sapphire substrate, consists of a central vacuum gap 
capacitor shunted by two coil inductors^’The bot¬ 
tom plate of the capacitor is split to create two cavity 
resonances, wi/27r = 8.89 GHz and a;2/27r = 9.93 GHz, 
named respectively the “measurement cavity” and “con¬ 
trol cavity”. The top plate of the capacitor is mechan¬ 
ically compliant, with a first harmonic mode of motion 
resonating at flm/27r = 14.98 MHz. Its motion tunes 
the resonance of both microwave cavities, with respec¬ 
tive vacuum optomechanical couplings gi/27r = 145 Hz 


and 52/271 = 170 Hz. Operated at a temperature of 
T = 30 mK, the equilibrium mechanical thermal occu¬ 
pancy is = 42 phonons and the mechanical relax¬ 
ation rate is T^/STr = 9.2 Hz. Both microwave cav¬ 
ities are strongly overcoupled to a single measurement 
port, setting their linewidths to Ki/27r = 1.7 MHz and 
K2/27r = 2.1 MHz. This coupling ensures that inter¬ 
nal dissipations contribute by less than 5% to the total 
linewidths, while maintaining a strongly resolved side¬ 
band regime, ^m /> 7- It also thermalizes the cav¬ 
ities to the shot-noise-limited input fields, maintaining 
throughout this work a thermal cavity occupancy well 
under our measurement noise floor, < O.I. 

We start by describing the QND measurement of the 
mechanical oscillator, cooled close to its ground state, in 
Fig.2. A cooling drive of strength T^ = 27r x 4.87 kHz = 
529 X Tm is applied at the lower mechanical sideband 
of the control cavity, = 022 — leading to a re¬ 
duced mechanical thermal occupancy rim- Simultane¬ 
ously, two drives of equal strength, T)" = T/", are applied 
close to the mechanical sidebands of the measurement 
cavity, acting back on the mechanical oscillator and in¬ 
creasing the total occupancy to = rim + n^a where 
riha = F)" /F^. Their frequencies can be optimally tuned 
to ujf' = ± Hm to perform a single mechanical quadra- 
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FIG. 2: Quantum nondemolition measurement, a, 

Measurement schematic. A cooling drive of strength = 
27r X 4.87 kHz = 529 x F^n is applied on the lower mechanical 
sideband of the control cavity = ui 2 — Hm). Two drives of 
equal strength F]" = F]^ are applied close to the mechanical 
sidebands of the measurement cavity. Their frequencies can 
be tuned to tu* = cui ± Qm to perform a single mechanical 
quadrature measurement (QND, in grey in b-d), or detuned 
by many mechanical linewidths away from that optimum to 
measure both mechanical quadratures (non-QND, in red in b- 
d). b, Mechanical noise spectrum (normalized, background 
subtracted), for FJ'/F^ = 0.9. c-d. Mechanical occupancy 
extracted from the measured spectra of the two-drives mea¬ 
surement (c) and cooling drive (d), for both the non-QND 
case (in red) and the QND case (in grey), as a function of the 
measurement strength T~[/T 2 . 


ture measurement (QND, in grey in Fig.2b-d), or detuned 
by many mechanical linewidths away from that optimum 
to measure both mechanical quadratures (non-QND, in 
red in Fig.2b-d). By monitoring the driven responses of 
both cavities^d9.20 tune very precisely the strength of 
each drive and measure all the mode frequencies and de¬ 
cay rates [Supp. Inf.]. We then acquire the noise power 
emitted by both cavities. In Fig.2b, we fix the measure¬ 
ment rate to F]" = F]^ = 0.9 xF^, and show the measured 
spectra, normalized to mechanical units [Supp. Inf.]. 

In the non-QND case, each drive measures both me¬ 
chanical quadratures, and the noise power of the thermo¬ 
mechanical sidebands are proportional to and + l 
for the anti-Stokes and Stokes scattering, respectively^^. 
Note that this sideband asymmetry^^’^^ provides a pri¬ 
mary calibration of the y-axis, in good agreement with 
the independently measured coupling strengths gi and 


g 2 - From each sideband, we extract the same mechani¬ 
cal occupancy shown as a function of measurement 
strength in Fig.2c-d. The measured quantum backac- 
tion scales ideally with the measurement strength, and 
we can extrapolate a mechanical thermal occupancy of 
= 0.15 ±0.05. 

We now tune the frequency of the measurement drives 
to the QND case. As shown in Fig.2b, the noise sideband 
of the cooling drive is unchanged. Indeed, that drive 
still measures both mechanical quadratures, accessing 
the same total mechanical occupancy as in the non-QND 
measurement (see Fig.2d). To the contrary, on the mea¬ 
surement cavity, the mechanical sidebands of each drive 
interfere with each other when brought into the QND 
case, leaving a single Lorentzian noise peak proportional 
to the variance of a single mechanical quadrature given 
by {Xi) = 2nm±l. The backaction has been evaded and 
placed on the orthogonal quadrature, conserving the total 
mechanical occupancy (see Eq.l and Eq.2). As expected 
the measured mechanical thermal occupancy is constant 
as a function of the measurement strength, and quantita¬ 
tively agrees with the occupancy inferred in the non-QND 
case. At the measurement strength of r]"/r^ = 2.44 we 
measure an evasion of the quantum measurement back- 
action by more than 13 dB. This demonstrates a QND 
measurement of a single mechanical quadrature at a rate 
much faster than the mechanical decoherence rate. 

This measurement scheme allows us to perform the to¬ 
mography of the mechanical state, described in Fig.3. In¬ 
deed, we can control the generalized mechanical quadra¬ 
ture being measured, , by simply rotating the relative 
phase between the measurement drives. As we expect 
the state to be Gaussian, the measurement of the sec¬ 
ond moment of the noise is sufficient to reconstruct its 
tomogram. A mechanical state prepared by simple side¬ 
band cooling is expected to have equal variances for each 
quadrature. In Fig.3c, we show data for an increased 
cooling strength F^ = 27r x 15.11 kHz = 1643 x F^ and 
a measurement strength r]"/r^ = 0.48. As expected the 
results of the QND measurements are phase independent. 
We measure a mechanical occupancy Um <0.1, demon¬ 
strating the QND measurement of a highly pure Gaussian 
state. 

We now apply this same tomographic measurement 
to verify the preparation of a non-classical state of mo¬ 
tion. We prepare a squeezed state by adding a drive 
of strength Fj at the upper mechanical sideband of the 
control cavity (uj^ = UJ 2 + D^). Again, as a squeezed 
state is still Gaussian, we simply measure the variance of 
the mechanical quadrature as a function of the mea¬ 
surement phase to fully characterize the state, as shown 
in Fig.3c for a squeezing strength FJ/F^ = 0.07. We 
resolve a minimum quadrature variance below vacuum, 
{Xi) = 0.78 ± 0.08. The spectra corresponding to the 
measurement of the squeezed and anti-squeezed quadra¬ 
tures are shown in Fig.3b. Finally, in Fig.3.d, we plot 
the variance of the squeezed and anti-squeezed quadra¬ 
tures, {Xi) and {X 2 ) respectively, as a function of the 
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FIG. 3: Tomography of a mechanical squeezed state, a, Measurement schematic. A pair of drives at tuj = 0^2 ± fim 
cool the mechanical mode to a squeezed bath, while a pair of drives at = oti ± Q,m measure the generalized mechanical 
quadrature given by their relative and tunable phase "F, allowing for the tomographic measurement of the squeezed state, 
b, Normalized noise spectra of the QND measurement (background subtracted), for the squeezed and anti-squeezed quadratures 
(r^/rj = 0.07) in blue and yellow respectively, compared to a spectrum measured without squeezing (F^ = 0) in grey, c. 
Measured quadrature variances as a function of phase for Fj = 0 in grey and F^/Fj = 0.07 in green. The spectra in b 
correspond to the blue and yellow dots, d. Squeezed and anti-squeezed quadrature variances, respectively in blue and yellow, 
as a function of the squeezing strength F^/Fj = 0.07. The solid lines are the theoretical prediction with no free parameters. 
The inset is a zoom-in for the low squeezing strength data, on a linear scale, where we can place the variance measured without 
squeezing (grey dot). In all figures, the black dotted line is the vacuum limit. 


squeezing strength. The solid lines in Fig.3c,d are the¬ 
oretical predictions from Eq.l and 2 without free pa¬ 
rameters, in reasonable agreement with the data at low 
squeezing strength. To further understand the deviation 
between data and theory at high squeezing strength, fu¬ 
ture experiments will investigate the effects of frequency 
noise of the cavities and the mechanical oscillator, resid¬ 
ual thermal cavity occupancy or phase instability of the 
drives. We emphasize that acquiring all the moments of 
the noise emitted by the measurement would allow us to 
reconstruct an arbitrary quantum state without making 
a Gaussian assumption^"^’^®. 

Looking forward, the introduction of stronger non- 


linearities, combined with reservoir engineering, would 
enable the preparation of more complex quantum 
states^*^, further motivating the use of mechani¬ 
cal systems as ultra-sensitive detectors and quantum 
memories®. Additionally, the inherent non-linearity 
of optomechanical cavities can act as a nearly-ideal 
mixing element, opening routes for innovative types 
of amplification^®, frequency conversion^^, and non¬ 
reciprocal behavior^®. 

Note: While preparing the manuscript we became 
aware of other work using a similar method for the me¬ 
chanical state preparation^®’®®. 
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Mechanical spectra and error estimation 

The mechanical state is extracted from the total microwave power spectral density, measured at room temperature 
with a spectrum analyser, following: 


PSD^ [g;] ^ _ / AnT ^ 

Pd [K2(l-2r;)2 + 4A2][K2+4(A±n„)2)] I^r2 + 4J2 + Y J xlp ^ ’ 

where PSD^ is the measured power spectral density at the upper or lower sideband of the drive in units of [W/Hz], 
Pd is the measured pump power in units of [W], A = — Wc is the drive detuning, 5 = w — ± is the detuning 

around the mechanical sideband, tOc is the microwave cavity frequency, Qm is the mechanical oscillator frequency, uid 
is the drive frequency, T is the total mechanical linewidth, k = Hint + Kext is the total cavity linewidth, rj = Kextjn 
is the coupling factor, and is the optomechanical coupling strength. Here, Uirnp is the measurement noise floor 
expressed in unit of mechanical quanta, and n = for the anti-Stokes scattering or n = + 1 for the Stokes 

scattering. 

From this measurement one can extract the mechanical displacement spectral density, Sxx, normalized by the 
variance of zero point fluctuation = h/{2mflm), where m is the mechanical oscillator’s mass. In Fig.2.b and 
Fig.S.b, we show the quantity S^xlx^p x F/4, allowing one to read the value of n off of the maximum value of the 
Lorentzian peak. We remove the noise floor B for readability, as the noise floor differs slightly between the control 
and measurement cavity. The transduction factor A and the linewidth F are measured independently, as explained 
in the next sections. 

In Figs.2 and 3, the error bar on the data correspond to the 90% confidence intervals, obtained from standard error 
propagation. The main source of uncertainties in our experiment are due, in order of greatest importance, to the error 
in the Lorentzian fits of the spectra, to the uncertainty in the estimation of gi ^2 and of the cavities coupling factor 

Calibrations and experimental details 

All the drives are filtered at room temperature and attenuated in the cryostat, ensuring that they are devoid of 
excess noise at the cavities frequencies^. 

The intrinsic mechanical relaxation rate F^ is measured from the mechanical ring-down time t^. In Fig.SI, we 
show the measured ring-down time as a function of the strength of a cooling drive on the control cavity, at a frequency 
uj^ = UJ 2 — flm, scattering light and damping the mechanical oscillator at the rate F^. We measure an intrinsic ring- 
down time Tm = 17.3 ms, corresponding to a relaxation rate F^ = 1/Tm = 27r x 9.2 Hz. This measurement also allows 
us to calibrate the cooling drive strength F^. 

The optomechanical coupling strengths for each cavity, 51,27 are calibrated from a fridge temperature sweep like in 
Teufel, et aP. The measured values are confirmed by the measurement of the sideband asymmetry in Fig.2. 

Finally the cavity frequencies and linewidths, the mechanical frequency, and every pump scattering rate, are 
independently measured from the driven responses (see next section). 

Scattering Parameters 

Before measuring the spectra shown in Fig.2 and Fig.3, we measure the driven responses of the system in presence 
of all the pump drives. They contain all the information about our system, except for the bath temperatures, and 
represent an important calibration. With up to four drives, two cavity modes and a mechanical mode, and with 
correction coming from finite sideband-resolution, the analytic solutions for the driven responses can be cumbersome, 
if not intractable. However, they can easily be solved numerically from the linear coupled equations of motion, 
following methods like in Andrews, et aP or Ranzani, et aV^. 

In Fig.S2 and Fig.S3 we show the driven responses corresponding respectively to the data shown in Fig.2 and 
Fig.3 in the main text. From the fit to these driven responses we extract the cavity frequencies and linewidths, the 
mechanical frequency, and every pump scattering rate. 


^ Lecocq F., et al, Resolving the vacuum fluctuations of an optomechanical system using an artificial atom, Nat. Phys. 11, 635 
(2015) 
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Cooling strength, /r„ 

FIG. SI: Mechanical ring-down time. We measure the mechanical ring down time Tm as a function of the strength of a 
cooling drive on the control cavity, at a frequency = uji — flm, scattering light and damping the mechanical oscillator at 
the rate F^. The solid line is the theoretical prediction from the total damping, Ftot = F™ + F^. We measure an intrinsic 
mechanical ring-down time of Tm = 17.3 ms, corresponding to a relaxation rate F^ = 1/rm = Stt x 9.2 Hz. 


^ Teufel J. D., et al, Sideband cooling of micromechanical motion to the quantum ground state, Nature 475, 359 (2011) 

® Andrews R. W., et al, Bidirectional and efficient conversion between microwave and optical light, Nat. Phys. 10, 321 (2014) 
Ranzani L., et al, A geometric description of nonreciprocity in coupled two-mode systems. New Journal of Physics 17, 023024 
(2015) 
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FIG. S2: Driven responses for the QND and non-QND measurement. Data corresponding to the measurement setup 
in Fig.2: the pump scattering rates are F^ = 27r x 4.87 kHz = 529 x F^ and Fj' = F^ = 0.9 x F^. a-b, Reflection coefficient, 
R, around each cavity resonance, c Transmission from the upper mechanical sideband of the red measurement drive to the 
lower mechanical sideband of the blue measurement drive, as a function of the input frequency. A signal sent at u} is received at 
uj — tjj” + uji . d Transmission from the upper mechanical sideband of the red cooling drive to the upper mechanical sideband 
of the red measurement drive, as a function of the input frequency. A signal sent at lj is received at w — + wj”. All the solid 

lines are from the same numerical simulation as a function of pump frequency, allowing for a very precise measurement of the 
cavity frequencies and linewidths, the mechanical frequency, and every pump scattering rate. 
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FIG. S3: Driven responses for the tomography of a mechanical squeezed state. Data corresponding to the measure¬ 
ment setup in Fig.3: the pump scattering rates are F^ = 2?! x 15.11 kHz = 1643 x Fm, 2 ~ and F^/F^ = 0.07. 
a-b, reflection coefficient, R, around each cavity resonance, c Transmission from the upper mechanical sideband of the red 
measurement drive to the lower mechanical sideband of the blue measurement drive, as a function of the input frequency. A 
signal sent at uj is received at ut — wj” d Transmission from the upper mechanical sideband of the red cooling drive to the 

upper mechanical sideband of the red measurement drive, as a function of the input frequency. A signal sent at w is received 
at a; — -|- u>^ . All the solid lines are from the same numerical simulation as a function of pump phases, allowing for a very 

precise measurement of the cavity frequencies and linewidths, the mechanical frequency, and every pump scattering rate. 




































